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Abstract. We consider a perturbed KdV equation: 

u + Uxxx - Quux = ef{x,u{-)), xeT, udx^Q. 

Jt 

For any periodic function u{x), let I{u) — /2(m), • • •) € be the vector, 

formed by the KdV integrals of motion, calculated for the potential u{x). Assuming 
that the perturbation ef{x,u{-)) is a smoothing mapping (e.g. it is a smooth function 
ef{x), independent from u), and that solutions of the perturbed equation satisfy some 
mild a-priori assumptions, we prove that for solutions u(t,x) with typical initial data 
and for ^ i < e^^, the vector I{u{t)) may be well approximated by a solution of the 
averaged equation. 
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We consider a perturbed Korteweg-de Vries (KdV) equation with zero mean-value 
periodic boundary condition: 

u + Uxxx — ^uux = ef{x,u{-)), X G T = R/Z, / u{x,t)dx = (). (0.1) 

Here ef{x,u{-)) is a nonlinear perturbation, specified below. For any p G M we denote 
by the Sobolev space of order p, formed by real- valued periodic functions with zero 
mean- value, provided with the homogeneous norm 1 1 ■ | |p. Particularly, if p G N we have 



HP = \ue L'^{T) : ||m||p < oo, / udx = o|, = / 



dPu 2 , 
dx. 



dxP 



nim,u2) = -(^{-^y\i,u2), (0.2) 



For any p, the operator ^ defines a linear isomorphism: ^ : — )• H^^^. Denoting 
by (^)~^ its inverse, we provide the spaces H^, p ^ 0, with a symplectic structure by 
means of the 2- form Vt: 

d_ 

'dx' 

where (■, ■) is the scalar product in L^(T). Then in any space H^, p ^ 1, the KdV 
equation flO.ip ^-n may be written as a Hamiltonian system with the Hamiltonian "H, 
given by 'H(m) = (^^ul + u^^dx. That is, KdV may be written as 

,. = |v«(,.). 

It is well-known that KdV is integrable. It means that the function space admits 
analytic symplectic coordinates v = (vi,V2,---) = \E'(m(-)), where = {vj,v_j) G M^, 
such that the quantities Ij = ||vjp, j ^ 1, are actions (integrals of motion), while 
ifj = Arg Vj, j ^ 1, are angles. In the (/, y?) -variables, KdV takes the integrable form 

/ = 0, ^ = W{I), (0.3) 

where W{I) G M.°° is the frequency vector (see [HEj). The integrating transformation 
\E', called the nonlinear Fourier transform, for any p ^ defines an analytic isomorphism 
m -.HP ^hP, where 



HP 



{v = (vi, V2, • ■ ■) : \v\l = 5^(27rjf P+V,f < oo, v, G M^ j G n}. 

It is well established that for a perturbed integrable finite-dimensional system, 

j = €/(/,¥.), ^ = W{I) + tg{I,^), e«l, 

where / G M", if G T*^, on time intervals of order the actions I{t) may be well 
approximated by solutions of the averaged equation: 

j = e(/)(J), (/)(J)= / /(J,^)rf^, 

provided that the initial data (/(O), (y9(0)) are typical (see [3l HJ [5l [6|). This assertion 
is known as the averaging principle. But in the infinite dimensional case, there is no 
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u — eu^ 



similar general result. In [TJ Ej, S. Kuksin and A. Piatniski proved that the averaging 
principle holds for the randomly perturbed KdV equation of the form: 

+ ,.„-6».. = v^„(M). /w. = /,,dx = 0, (0.4) 

where the force rj is a white noise in t, is smooth in x and is non-degenerate. Our 
goal in this work is to justify the averaging principle for the KdV equation with 
deterministic perturbations, using the Anosov scheme (see [3]), exploited earlier in the 
finite dimensional situation. The main technical difficulty to achieve this goal comes 
from the fact that to perform the scheme one has to use a measure in the function 
space which is quasi-invariant under the flow of the perturbed equation (it is needed 
to guarantee that a small 'bad' set which we have to prohibit for a solution of the 
perturbed equation at a time t > corresponds to a small set of initial data). For a 
reason, explained in Section 3, to construct such a quasi- invariant measure we have to 
assume that the perturbation ef is smoothing. More precisely, we assume that: 

Assumption A. (i) For any p ^ 0, the mapping deflned by the perturbation in (10. ip : 

V -.HP ^ HP+^'\ u ^ /(x, m(-)), (0.5) 

is analytic. Here Co > 1 is a constant. 

(ii) For any J9 ^ 3 and T > 0, the perturbed KdV equation (10. ip with initial data 

m(0) = Mo G i/P, 

has a unique solution u(t,x) G in the time interval [— Te^^Te^^], and 
\\u{t)\\p^C{p,\\uo\\p,T), \t\ ^Te-\ 

We are mainly concerned with the behavior of the actions I{u{t)) G M^jf for \t\ < e^^. 
For this end, it is convenient to pass to the slow time t = et and write the perturbed 
KdV equation (10. ip in the action-angle coordinates {I,(p): 

^ = F(J,(^), ^=e-'WiI) + GiI,^). (0.6) 

Here / G M°°, G T°° and T°° := {6 = (0,)i^i, Oi G T} is the inflnite-dimensional torus, 
endowed with the Tikhonov topology. The two functions F{I, ip) and G{I, ip) are the 
perturbation term e/, written in action-angle variables, see below (II. 3p and (II. 4p . The 
corresponding averaged equation is 

{F)iJ), {F)iJ)= [ FiJ,ip)dip, (0.7) 



dr 

where dip is the Haar measure on T°°. It turns out that the (10.71) is a Lipschitz equation, 
see below (I4.17p . We denote by the image of the space under the action- mapping 



1. ,2 



Clearly, / = 7ij{v) G C h^j, where h^j is the weighted /^-space 

oo 

G M~ : |/|,. = |/|, = 2 J](27rjf < oo}. 
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and h^j_^ is its positive octant, h^j_^ = {I E h^j : Ij ^ 0, Vj'}. This is a closed subset of h^j. 

For any 9 = (6'j)j^i G T°°, let us denote by $51 the linear operator on the space of 
sequences (vi, V2, • ■ ■) G /i^ which rotates each component Vj G by the angle 9j. 

Definition 0.1 A Gaussian measure fi on the Hilbert space is said to be (0 -admissible 
(where Co > 1 is the same as in assumption A), if the following conditions are fulfilled: 

(i) It is non-degenerate and has zero mean value. 

(ii) It has a diagonal correlation operator (vi,V2,---) i-> (criVi, (T2V2, ■ ■ ■), where 
every aj > 0, Xlj^ji '^j < ^ = 0(1). In particular, /i is invariant under the 
rotations $51. 

Such measures can be written as: 

j=i J 3 

where d\j, j ^ 1, is the Lebesgue measure on (see [HIIU]). Clearly, they are invariant 
under the KdV flow fl03|) . 

The main result of this work is the following theorem: 

Theorem 0.2. Fix any p ^ 3 and T > 0. Let the curve u'(t) G H^, \t\ ^ e^^T be a 
solution of equation f lO.ip and f "(r) = \I^(M''(e"^r)), r = et, |r| ^ T. If assumption A is 
fulfilled and /i is a (^q- admissible Gaussian measure on h'^, then 

(i) For any p > 0, there exists a Borel subset of /i^ and ep > such that 
lim^^o yu(^^ \ rp = 0, and for e ^ we have 

\Iiv'iT))-Jir)\,^p, for IrKf, t;^(0) G F^, (0.9) 

where J{t), \t\ ^ T, is a solution of the averaged equation (10. 7p with the inital data 
J(0) = njivm)- 

(ii) There is a full measure subset F^ of with the following property: 

If t'^(O) G F^, then for any ^Ti <T2 the image fij,^ j,^ of the probability measure 
(T2 — Ti)^^dT on [Ti,T2] under the mapping r h-)- (y9(ti''(r)) G T°° converges weakly, as 
e — )■ 0, to the Haar measure dip on T°°. 

The assertion (ii) of the theorem means that for any bounded continuous function 
g{^) on T~, 

7fr~r [_ g{v{v'{T)))dT [ g{^)d^, e 0. 

-'2 ^ 1 JTi JT°= 

In particular, we have 

Proposition 0.3. The assumption A holds if in (10. ip / = f{x) is a smooth function, 
independent from u. 

It is unknown for us that if the result of Theorem 0.2 remains true for equation 
(10. ip with non-smoothing perturbations, e.g. if the right hand side of equation (10. ip is 
eUxx or —eu. So we do not know whether a suitable analogy of the result in [71 18] holds 
true if in equation (10. 4p the noise rj vanishes. 
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The paper has the following structure: Section 1 is about the transformation which 
integrates the KdV and its Birkhoff normal form. In Section 2 we discuss the averaged 
equation. We prove that the (^o-admissible Gaussian measures are quasi-invariant under 
the flow of equation (10. ip in Section 3. Finally in Section 4 and Section 5 we establish 
the main theorem and Proposition 0.3. 

Agreements. Analyticity of maps Bi — t- B2 between Banach spaces Bi and B2, which 
are the real parts of complex spaces B^ and B2, is understood in the sense of Frechet. 
All analytic maps that we consider possess the following additional property: for any R, 
a map extends to a bounded analytical mapping in a complex > 0)-neighborhood 
of the ball {\u\b^ < R} in B^. 

Notation. We use capital letters C or C{ai,a2, ■ ■ ■) to denote positive constants that 
depend on the parameters ai, 02, ... but not on the unknown function u. We denote 
Z^o = G Z, n ^ 0}. For an infinite-dimensional vector w = {wi,W2, ■ ■ ■) and any 
n G N we denote = {wi, . . . , w„, 0, 0, . . .). We often identify with a corresponding 
n-vector. 



1. Preliminaries on the KdV equation 

In this section we discuss integrability of the KdV equation (0.1)e=o- 



1.1. Nonlinear Fourier transform for KdV 

We provide the L^-space with the Hilbert basis {e^, s G Z \ {0}}, 

\/2 cos(27rsx) s > 0, 
-\/2 sin(27rsx) s < 0. 

Theorem 1.1. There exists an analytic diffeomorphism \1' : and an analytic 

functional K on of the form K{v) = K{I{v)), where the function K{I) is analytic in 
a suitable neighborhood of the octant h^j^ in h^j, with the following properties: 

(i) The mapping defines an analytic diffeomorphism \1/ : i— )■ h^, for any 
p G Zj>o- This is a symplectomorphism of the spaces {Hp,Q) (see (10.21) and {hP,uj2), 
where 002 = ^1 ^'^k A dv-k- 

(ii) The differential d^{0) takes the form '^UgCg ^ v,Vs = \2tts\^^^'^Us. 

(iii) A curve u G C^(0, T; is a solution of the KdV equation (lO.ip .-n if and only 
if v(t) = \l/(-u(t)) satisfies the equation 

[^1 ~o)w^^^^^'' V, = (t;„t;_,■)GM^ JGN. (1.1) 
Since the maps and are analytic, then for m = 0, 1, 2 . . ., we have 

\\d^^{u)\U < PmiMm), \\d^^-\v)\\^ ^ Qmi\v\m), J = 0, 1,2, 



where Pm and Qm are continuous functions (cf. the agreements). 
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We denote 

iy(/) = (1^1,1^2,...), = |^(/), A; = 1,2,.... 

Lemma 1.2. For any n G N, if In+i = In+2 = ■ ■ ■ = 0, then 

, f ,dWi. \ , 
det(^(^)i^i,K„j t 0. 

Let 1°°^ be the Banach space of all real sequences / = (/i, • • •) with the norm 
= supr2^"'"|/„| < oo. 

Denote k = (k„)„^i, where k„ = (2?™)^. 

Lemma 1.3. The normalized frequency map 

W -.1^ W{I) = W{I) - K 

is real analytic as a map from h]_^ to I'^i- 

The coordinates v = "^{u) are called the Birkhoff coordinates, and the form (1.1) of 
KdV is its Birkhoff' normal form. See [1] for Theorem 1.1 and Lemma 1.3. A detailed 
proof of Lemma 1.2 can be found in [2]. 

1.2. Equation Ii0.1\) in the Birkhoff coordinates. 
For /c = 1, 2 ... we denote: 

where = = (vi, V2, . . .). Let u(t) be a solution of equation (10. ip . We get 

V, = d^kiu)iefix,u) + Viu)), k^l, (1.2) 

where V(u) = —u^xx + Guu^. Since Ik{v) = ll^&fcp is an integral of motion of KdV 
equation (lO.ip .-n. we have 

4 = e{d^k{u)f{x, u), Vfc) := eFk{v). (1.3) 

Here and below (■, ■) indicates the scalar product in R^. 

For k ^ 1 define (pk = arctan(— ) if ^ 0, and y^fc = if = 0. Using equation 
(1.1), we get 

^k = WkiI) + e\Yk\-\d-^k{u)fix,u),vi), if Vfc^O, (1.4) 

where = {—V-k,Vk)- Denoting for brevity, the vector field in equation (1.4) by 
Wk{I) + eGk{v), we rewrite the equation for the pair (7^, ipk)ik ^ 1) as 

ik{t) = eFk{v) = eFk{I,^), 
ipkit) = Wk{I) + eGkiv). 

We set 

F(J,^) = (Fi(/,y.),F2(J,y.),...). 
In the following lemma Pk and P^ are some fixed continuous functions. 
Lemma 1.4. For k,j G N, we have for any p ^ 
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(i) The function Fk{v) is analytic in each space h^. 

(ii) For any p ^ 0, 5 > 0, the function Gkiv)x{ik^s} is bounded by S~^^'^Pk{\v\p). 

(iii) For any 6 > 0, the function ^^{I,ip)x{ij^5} is bounded by 6~^^^Pi{\v\p). 

(iv) The function |^(/, ^) is bounded by P^(|t>|p), and for any n eN and (Ji, . . . , /„) G 
MJ^, the function F^^Ii, ipi, . . . ,In, ipn, 0, . . .) is analytic on T". 

Proof. Items (i) and (ii) follow directly from Theorem 1.1. Items (iii) and (iv) follow 
from item (i) and the chain-rule: 

dFk rr^^dFk , . OF, 



- ^''H9^^°^(^^■)-9^^^"(^- 



dij y J \ Qy. 

From this lemma we know that equation (11. 5p may have singularities at dh^j_^. We 
denote 

Uj:hP^ hP, Ui{v) = I{v), 
n,,^ : /iP -> /,P X T°°, Uj^^iv) = (liv), ^{v)). 
Abusing notation, we will identify v with (/, v?) = ^i,ip{v)- 

Definition 1.5. For p ^ 3, we say that a curve {I{t),(f{t)), \t\ ^ T, is a regular 
solution of equation (11.51) . if there exists a solution u(t) E H'^ of equation (10.11) such 
that u{t) E HP and 

n,,^(vl/(M(t))) = (/(t),v9(t)), \t\<:T. 

If (/(t), f(t)) is a regular solution of (II. 5p and |/(0)|p ^ Mq, then by assumption A 
we have 

\nt)\p = mil ^ Cip,Mo,T), \t\ ^ Te~\ (1.6) 
2. Averaged equation 

For a function / on a Hilbert space if, we write / E Lipioc{H) if 

\fiu,)- fiu2)\^PiR)\\u,-U2\l \\Ui\\,\\u2\\^R, (2.1) 

for a suitable continuous function P which depends on /. Clearly, the set of functions 
Lipioc{H) is an algebra. By the Cauchy inequality, any analytic function on H belongs 
to Lipioc{H) (see agreements). In particularly, for any k ^ 1, 

Wk{I) E LipiocifiPj), p^l, and Fk{v) E Lipioc{hP), p ^ 0. 

In the further analysis, we systematically use the fact that the functional F^iv) only 
weakly depends on the tail of the vector v. Now we state the corresponding results. Let 
/ E LipiocihP) and v E h^^ , pi > p. Denoting by 11*^, M ^ 1 the projection 

: /,o ^ /^o^ (^^^ . . .) (vi, . . . , VA/, 0, . . .), 
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we have \v — Yi^^v\p ^ {2txM)^'^'^^~p^\v\p^. Accordingly, 

\f{v) - /(n*Ml ^ P{\vU{27:M)-^P^-^\ (2.2) 

The torus T^^ acts on the space T^Mh^ by hnear transformations ^Omi ^ T^'^ , 
where : {Im,Vm) ^ {Imi^m + Om)- Similarly, the torus T°° acts on /i° by linear 
transformations $e : (/, Lp) i— )■ + 9) with 9 E T°°. 

For a function / G Lipioc{hP) and a positive integer we define the average of / 
in the first angles as the function 

{f)N{v)= [ f({^e^®ld)iv))d9^, 

and define the averaging in all angles as 



{f)iv) = / fiMv))de, 

J fee 

where d9 is the Haar measure on T°°. The estimate fl2.2p readily implies that 

Let V = (/, (p), then (/) jy is a function independent of ■ ■ ■ , v^at, and (/) is independent 
of ip. Thus (/) can be written as (/)(/)• 

Lemma 2.1. (See [7]). Let / G LipiodhP), then 

(i) The functions {f)N{v) and (/)(f) satisfy (12. ip with the same function P as f and 
take the same value at the origin. 

(ii) These two functions are smooth (analytic) if / is. If / is smooth, then (/)(/) is a 
smooth function with respect to vector (Ji, . . . , Im), for any M. If /(t>) is analytic 
in the space h''', then (/)(/) is analytic in the space h^j. 

We recall that a vector w G M" is non-resonant if 
oo-k^O, VA; G \ {0}. 
Denote by C'^^^(T") the set of all Lipschitz functions on 

Lemma 2.2. Let / G C*^^^(T") for some n E N. Then for any non-resonant vector 
G M" we have 



1 r 

hm - / /(xo + u;t)dt = (/), 



uniformly in G T". The rate of convergence depends on n, u and /. 

Proof. Let us write f{x) as the Fourier series f{x) = ^ fk^^^^ ■ Since the Fourier series 

of a Lipschitz function converges uniformly (see [H]), for any e > we may find R = 



such that 



^ I for all X. Now it is enough to show that 
^ fR{xo + cot)dt-fo\^^, VT^T,, (2.3) 
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for a suitable T^, where = J2\k\<Rfk^^''^- Observing that 



1 







T\k ■ uj\ 

for each nonzero k. Therefore the l.h.s of (2.3) is smaller than 

2 / . ^-l 



— ( inf I A; ■ cjI ) \fk\. 

\k\s:R 

The assertion of the lemma follows. □ 



3. Quasi-invariance of Gaussian measures 

Fix any integer p ^ 3, and let be a Co-admissible Gaussian measure on the Hilbert 
space hP. In this section we will discuss how this measure evolves under the flow of 
the perturbed KdV equation (10.11) . We follow a classical procedure based on finite 
dimensional approximations (see e.g. [T2l flO] ). 

We suppose the assumption A holds. Let us write the equation (10. ip in the Birkhoff 
normal form, using the slow time r = et: 

Av^. = e-i:ri^,(/)v, + X,(t;), jGN, (3.1) 



where = {Xj, X^jY G and J 



-1 


For any n G N, we consider the 2n-dimensional subspace TCnih^) of with 
coordinates = (vi, . . . , v„, 0, . . .). On Hn{h^), we define the following finite- 
dimensional systems: 

= e-^JWj{I{u''))uj + X,(w"), 1 ^ J ^ n, (3.2) 

dr 

where Uj = {uj, (jJ-jY ^ and = (cJi, . . . , cD^, 0, . . .) G TTn{hP). 

We denote X"(f") = (Xi(t;"), . . . , X„(t;"), 0, . . .) and X{v) = (Xi{v) , . . .) . By 
assumption A and Theorem 1.1, for any p ^ the mapping 

X:hP ^ W^^'^ V ^ X{v) is analytic. (3.3) 

Theorem 3.2. For any T > 0, uj"'{-) converges to v{-) as n — t- oo in C([— T, T]; /i^), 
where f(-) and w"(-) are, respectively, solutions of (13. ip and ( 13. 2p with initial data 
v{0) G hP and w"(0) = f'^(O) G 7r„(/iP). 

Proof. Fix any Mq > 0. From (II. 6p we know that there exists a constant Mi such 
that if \v{0)\p ^ Mo, then 

|t;(r)|p^Mi, rG[0,r]. (3.4) 

The equation (13. 2p yields that 

-\ioX = 2^ji+2^a;, ■ X,(a;") := x"(u;"). (3.5) 
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We define 

oo 

By (13.31) . we know that there exists a constant Ci > such that 

Ix'^MI ^ Ci, |w"|p^2Mi, VneN. (3.6) 
Denote f = Mi/Ci, then if |u;"(0)|p ^ Mq, then 

|a;"(r)|p ^ 2Afi, rG[-f,f], G N. (3.7) 
Lemma 3.3. In the space C([— f, r], h^^^), we have the convergence 

ci;"(-) —7- f (■) as n — 7- oo. 

Proof. Denote C,j = — cDj, = /(f) and J^^n = /(w"). Since ^/vj = v^, using 
equations (13. ip and (13. 2p . for 1 ^ j ^ n, we get 

^iSr = 2(0)*[e-V(WG-(/.)v, - iy,(/..)cl;,) + X,(f ) - X,(a;")] 

= 2e-i[l^,(4) - W,-(/.n)]v, ■ (a;,)^ + 2(g)* ■ (X,-(f ) - X,(a;")). 
By Lemma 1.3 and Cauchy's inequahty, we know that 

W,{I{v)) - l^,(/(u;"))| ^ C2{Mi)j\v-uj\.i. 

Using (13. 3 p we get that 

^\v- uX-i ^ C,{e, Ml) \v - uX-i + an{v), T e i-f, f], 

where 

oo 

j=n+l 

Obviously, a„(f) — )■ as n — )• oo uniformly for \v\p ^ Mi. 

The lemma now follows directly from Gronwall's Lemma. □ 

Lemma 3.4. If ^"(0) — f (0) strongly in and — i- r, r„ G [— r, f], as n — oo, then 
lim \v{t) - a;"(r„)|p = 0. 

n— >oo 

Proof From (13.51) we know that for any r„ G [— r, f]. 



X"(a;"(s))ds. 



Since a;"(0) — )■ f (0) strongly in h^, then using (13. 3p and Lemma 3.3 we get 
\v{t)\1 ^ hminf |u;"(r„)P ^ limsup |a;"(r„)P 



n— >oo 

Tn 



limsup(K(0)|J+ / x"(c^"(5)M5) = K0)i;+ / X(^(^))rfs 
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Therefore, lim„_j.oo \^^ijn)\p = |'y(''")|p- Since a;"(r„) — )■ v{t) in the space /i^"^ as n — )■ oo, 
then the required convergence follows. □ 

Lemma 3.5. In the space C([— f , f], /i^), Co'"(-) — )■ f (■) as n — )■ oo. 

Proof. Suppose this statement is invalid. Then there exists 5 > and a sequence 

{r'^}„gN C [— r, f ] such that 

|u;"(r") -t;(r")|p ^ 5. 

Let {r'^'^jfegN be a subsequence of the sequence {r"}„gN converging to some r° G [— r, f]. 
But t'(r"'=) —7- v{t^) in as k oo, and using Lemma 3.4, we can get a;"'=(r"'=) — )■ f (r") 
as — 7- oo in h^. So we get a contradiction, and Lemma 3.5 is proved. □ 

If T ^ f, the theorem is proved, otherwise we iterate the above procedure. This 
finishes the proof of Theorem 3.2. □ 

Let denote the flow determined by equations (13. ip in the space h^, and 

B;{M) := {vehP-. \v\p ^ M}. 

Theorem 3.6. For any Mq > and T > 0, there exists a constant C > which depends 
only on Mq and T, such that if A is a open subset of Bp{Mo), then for r G [0,T], we 
have 

Proof. From (11. 6p we know that there is constant Mi which only depends on Mq and 
T, such that if v{0) G B^{Mo), then 

v{r) G i?;(Mi), |r| ^ T. (3.8) 

For any G N, consider the measure /i„ = 7r„ o /i on the subspace 7r„(/iP). Since is a 
Co-admissible Gaussian measure, by (10.81) /i„ has the following density with respect to 
the Lebesgue measure: 

" 1 " ■l+2p| 12 

6.(^'^) := (27r)-"n(2vrjr+V'^^P{-2E ^ ^ >• 
i=i i=i ^ 

Let 5^ be the flow determined by equations (13.20 on subspace TTnih^)- For any open set 
An C 7r„(i?p(Mo)), due to Theorem A in the appendix, we have 

^f,n{S:{An)) 



n 



^/VjXj + V-jX_j dXj dX 



c"(t;")6„(t;")c/t;'^ 
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Since j~'^°/<Jj = 0(1), using fl3.3p and the Cauchy's inequality, there exists a constant 
C which depends only on Mi, such that 

|c"(t;")|^C, t;" G 7r„(5;(Mi), Vn G N. (3.9) 

We have 

e-^>„(A„) ^ finiS^iA^)) ^ e^>„(A„), (3.10) 

as long as 5;(A„) C 7r„(S;(Mi)). 

Since /i„ convergences weakly to fi, the theorem follows from (13. Sp . (I3.10p and 
Theorem 3.2 (see [121 [lO]). □ 

4. Proof of the main theorem 

In this section we prove Theorem 0.2 by developing a suitable infinite-dimensional 
version of the Anosov scheme (see [31 HJ [5l [6]), and by studying the behavior of the 
regular solutions of equation (11. 5p and the corresponding solutions of (10. ip . We fix 
p ^ 3. Assume m(0) = Mq G H^. So 

n,,^(vl/(Mo)) = (Jo, ^o) e h% X T~, 3. (4.1) 

4-1. Proof of the assertion (i) 
We denote 

Bl{M) = {/ G : ^ M}. 

Without loss of generality, we assume that T = 1 and t ^ 0. 
Fix any Mq > 0. Let 

(/o,<^o) G5p^(Mo) xT°° :=ro, 

that is, 

Vo = ^{uo) G i?;(v^). 

Let {I(t),ip{t)) be a regular solution of the system (11.50 with (/(0),v9(0)) = {Io,(Pq). 
Then by f ll.6p . there exists Mi ^ Mq such that 

/(t)G5;(Mi), tG[0,e-i]. (4.2) 

By the definition of the perturbation we know that 

\F{I,^)\,<:Cm„ V(/,<^)Gi?;(Mi)xT~, (4.3) 

where the constant Cmi depends only on Mi. 

We denote = (/i, . . . , 0, 0, . . .), (^"^ = ((^i, . . . , yp^, 0, 0, . . .), and iy"(J) = 
iW,{I), . . . , W„(/), 0, 0, . . .), for any m G N. 

Fix no G N. By (2.2), for any p > 0, there exists mo G N , depending only on uq 
and p, such that if m ^ mg, then 

|Ffc(J,^)-Ffc(/™,^"^)| ^p, V(/,^)Gi?;(Mi)xT°°, (4.4) 



(4.6) 
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where A; = 1, ■ ■ ■ , no- 

From now on, we always assume that 

(/,<^)g5;(Mi)xT~, i.e. veB;{^/M[). 

By Lemma 1.4, we have 

,4.5) 

From Lemma 1.3 and Lemma 2.1, we know that 

\W,{I) - W,{I)\ ^ Ci(j, Mi)|/ - 

|(Ffc)(/) - ^ Ci(A;,j,Mi)|/-/|i. 

By (2.1) we get 

|F,(/™°,y."^») -F,(J"^°,^"^»)| ^ C2{k,mo,M,)\v"'^-v"'% (4.7) 

where | ■ | is the maximum norm. 
We denote 

Cm^ = ^0 ■ max{Co, Ci, C2 : 1 ^ j ^ mo, 1 ^ ^ ^o}. 

Below we define a number of sets, depending on various parameters. All of them also 
depend on mo and uq, but this dependence is not indicated. For any 6 > 0, and To > 0, 
we define a subset E{6,Tq) C Bp{Mi) as the collection of all / G Bp{Mi) such that for 
every (f G T°° and any T ^ Tq, we have 

^ j\Fk{I"'\^"^^ + W"^%I)t) - {Fk){I'^°)]dt\^6, (4.8) 

for k = 1, ■ ■ ■ , uq. Let 5* be the fiow generated by regular solutions of the system (1.5). 
We define two more groups of sets. 

S{t) = S{t,e,6,To,I,^) := {ti G [0,t] : 5^(7, <^) ^ E{5,To) x T~}. 

N{f) = N{f, e, 6, To) := {(/, </^) G Fq : Mes[S(e-\ e, 6, To, /, y.)] ^ f }. 

Here and below Mes[-] stands for the Lebesgue measure in M. 

Clearly, E{6,Tq) is a closed subset of B^^Mi) and S(t,6,TQ, I,ip) is a open subset 
of [0,t]. The following result is the main lemma of this work: 

Lemma 4.1. For k = 1, . . . , hq, the /^-component of any regular solution of (11. 5p with 
initial data in N{T, e, 6, To) can be written as: 

Ik{t) = 4(0) + e [\Fk){I{s))ds + S(t), 
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where for any 7 G (0, 1) the function is bounded on [0, i] by 

4602"^' { [2(7 + 2ToCAf,e)i/2j (To + T + e'') 

+ [ y/2 + ^oCM,e + {^;^ + )^{Io + I+e )| 

+2eCM,f + 2p + 26 + 2eCM, (Tq + T). 

Proof For any (J, (^) e N{T), we consider the corresponding set It is composed 
of open intervals of total length less than min{T, t}. Thus at most [T/Tq] of them have 
length greater than or equal to Tq. We denote these long intervals by (oj, foj), 1 ^ i ^ d, 
d ^ T/Tq and denote by C(t) the complement of Ui^i^rf(aj, 6j) in [0,t]. 
By we have 

Fk{I{s),^{s))dt= [ Ffe(J"^«(s),<^"^°(s))ds + ei(t), 



'0 JC(t) 

where 16(^)1 < C^^T + pt. 

The set C{t) is composed of segments [6j_i,aj] (if necessary, we set 60 = 0, and 
ad+i = t). We proceed by dividing each segment [bi^i,ai] into shorter segments by 
points fj, where hi = t\ < t2 < ■ • • < t\. = ai. The points lie outside the set S{t) 
and To ^ t^^^^ — ^ 2Tq except for the terminal segment containing the end points Oj, 
which may be shorter than Tq. 

This partition is constructed as follows: 



If — ^ 2To, then we keep the whole segment with no subdivisions. {t\ = bi-i 

If ttj — > 2To, we divide the segment in the following way: 

a) If + 2Tq does not belong to S(t), we chose = ^j-i + 2To, and continue 
by subdividing [^2,01]; 

b) if bi-i + 2To belongs to S(t), then there are points in + To,6i_i + 2To] 
which do not, by definition of We set ^2 equal to one of these points and 
continue by subdividing [t2 5'^i]- 

We will adopt the notation: /i* = fj_^_^ — f- and j) = [t]-,t]j^-\\- So 

d rii — l 

= U U ^(''•^■)' ^0 ^ ^5 = l^^^'-^')! ^ 2To, J ^ - 2. 

2=1 i = l 

By its definition, C{t) contains at most [T/Tq] + 1 segments a-j], thus C{t) 
contains at most [T/Tq] + 1 terminal subsegments of length less than Tq. Since all 
other segments have length no less than Tq and t ^ the number of these segments 
is not greater than [eTo]~^. So the total number of subsegments s{i,j) is bounded by 
l + [(eTo)-i] + [f/ro]. 
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For each segment s{i,j) we define a subset A(i, j) of {1, 2, ■ ■ ■ , mo} in tlie following 

way: 

/gA(z,j) ^ 3tes{i,j), Ii{t)<-f. 
If / e A, then by (l43l) we have 

\Ii{t)\<2ToCM,e + ^, tes{i,j). (4.9) 
For / = (Ji, I2, . . .) and ip = {ipi, ip2, . . .) we set 

where = {ipi, 02, ■ ■ •) and / = (Ji, I2, ■ ■ ■) are defined by the following relation: 

If leA{i,j), then h = 0, 0i = 0, else h = h, 0i = ^pi- 

We also denote Xij{I,(p) = \ij{(p)) and when the segment s{i,j) is clearly 

indicated, we write for short Xij{I,ip) = {1,0)- 

Then on s{i,j), using (14.71) and (14.91) we obtain 



Fj /™" (s) , (s) - Fj A,,, (/"« (s) , y^-o (s) ) 



ds 



1/2 



ds 



/™"(s)-\,(/™»(s 

^2ToC;™(7 + 2ToC,,,e)i/2. 
In Proposition 1-5 below, k = 1, . . . ,no. 
Proposition 1. 

Ffc(/™«(s),(/."^«(s))rfs = 5^ [ Ffc(/-«(t;.),¥^"^''(s))rfs + 6(t) 



(4.10) 



C(t) 



where 



161 ^ 4C;;«;'"« [(7 + 2ToCM,tf" + 7-'/'ToCM,eJ {T, + f + 6"^). (4.11) 
Proof. We may write ^(t) as 



ds 



For each s{i,j), we have 

Ffc f /™« (s) , (^"o (s) ) - Ffc ( (t}) , <^-« (s) 



ds 



no, mo 



ds 



(4.12) 
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We replace the integrand Ffc(J'^o, by Ffc(J'^o, (^""o). Using fOUD and <KV]^ we 
obtain that 

The inequahty M.lip follows. □ 

On each subsegment s{i,j), we now consider the unperturbed linear dynamics ^p]{t) 
of the angles ip""" E T"*" : 

Proposition 2. 

5^ / FJl^"{t)),v"^^{s))ds = J2 [ FJl"^^{t)),v){s))ds + Ut), 



where 



1^3(^)1 ^ 4C;:?;™°(7 + 2ToCM,e)^/'(To + f + e-^) 



'Ml 



V 7 

Proof. For each s{i,j) we have 



+ 



(To + T + e 



-1^ 



(4.13) 



ds 



ds'ds 



{hi) Ji 



'-'Ml 



ds'ds 



'-'Ml 



7-'/'e(s-t;) + -CM,e(s-t;)^ 



^ ^Mi [ H ^ 



(is 



Here the first inequality comes from equation (II ■41) . and using (14. 5 p and (14. 6 p we can 
get the second inequality. The third one follows from (14. 3p . 
Using again (14. 5p . we get 



F, (/™o (t;. ) , (3) ) - F A,, (/-o [t] ),cpUs)) 



ds 



'-'Ml 



A,,(^-»(s)-^}(s)) 



ds 



Therefore (I4.13P holds for the same reason as (14. lip . □ 

We will now compare the integral Fk{I"^° (t^j) , ip'^j{s))ds with the average value 
(F,(J-o(t;.)))/,;.. 
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where 



\^,(t)\^- + 2CMATo + f). 



(4.14) 



Proof: We divide the set of segments s{i,j) into two subsets Ai and A2. Namely, 
s{i,j) G Ai if h^j ^ To and s{i,j) G A2 otherwise, 
(i) s{i,j) G Ai. In this case, by (14.81) . we have 



FJn^it)),^){s))-{F,)(l-^it)) 



ds 



^ 5h). 



So 



(ii) s{i,j) G A2. Now, using (14. 3 p we get 



F, ( (t;.) , ^}(s) ) ds - (Ffc) ( (t;.) ) /i;. 



s{i,j)eAi 



Since Card(A2) ^ (1 + T/Tq), then 



E 

s(ij)GA2 



F ( (tj ) , (s) ) rfs - (F,) ( (t;. ) ) ^ 2Cm, (T + Tc 



This imphes the inequahty (I4.14p . □ 
Proposition 4. 



Y.h){Fu){i'-'{ti) 



C7(t) 



{F,)(l^\s)]ds + Ut) 



where 



\Ut)\ ^ AtCM^CZr^oiT, + f + e-'). 
Proof Indeed, as 



(4.15) 



1^5(^)1 =|E 



ds 



jF,)(/-«(s))-(F,)(J-(t}))_ 
using (14. 3 p and (14.60 we get 

Finally, 
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Proposition 5. 

^^^^ (F,) (/™" {s)) ds = J\f,) (/(s)) ds + Ut) , 

and 1^6(^)1 is bounded by Cm^T + pt. □ 

Gathering the estimates in Propositions 1-5, we obtain 



= 4(0) + e f (Ffe)(/(s))ds + 5(t), 



where 



\m\<ej2m)\ 



1=1 



2(7 + 2roCM,6) V2 + + ToCM.e 



r 

^2, 



+ + ^^)] (To + T + e-) + 2eCM.T, 

+ 2p + 25 + 2eCMi(To + r), te[0,^]. 
Lemma 4.1 is proved. □ 

Corollary 4.2. For any p > 0, with a suitable choice of p, 7, 5, To, T, the function 
|S(i)| in Lemma 4.1 can be made smaller than p, if e is small enough. 
Proof. We choose 



with 



a „ 1 

1---(T>0, 0<(7<-. 

2 ' 2 



Then for e sufficiently small we have 

\E{t)\<p. □ 

On the Hilbert space h^, we adopt a Co-admissiblc Gaussian measure /i. Define 
corresponding measures /// = H/o// and /x/ ,^ = H/^^o// in the spaces h^_^_ and x T°°. 
Lemma 4.3. The measure /ii^^ is a product measure d/ij^^ — d/ijdip, where dip is the 
Haar measure on T°°. 

Proof: Since the measure p, is invariant under rotations ^0, the o dp is a measure 
on T°°, invariant under the rotations. So this is the Haar measure d(p. Gonsequently 
the image of the measure p/,^ under the natural projection (/, t-^ 93 is dip. Since 
the spaces and T°° are separable, then for (p e T°° there exists a Borel probability 
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measure 'n'^p{dl) on h^j^ such that fij^^ = 'n'^{dl)dip. That is, for any bounded continuous 
function /(/, (p), we have 




(see e.g. [9]). For any Q G we have 

f{I,^ + 9)7T^{dI)d^ = j j f{l,^)7i^^e{dl)d^. 

Integrating in d6 we see that 

l^i,^{,dldip) = dfi'{dl)dip, 

where dfi'{dl) = J^^c 'ne{,dl)dip. We must have dii' = dfij, and the assertion of the 
lemma is proved. □ 

The two lemmas below deal with the sets E and N, defined at the beginning of this 
section. 

Lemma 4.4. For any 6 > 0, limro^oo /i/(5p(Mi) \ E{6,To)) = 0. 
Proof: From the definition of E{6,To), we know that 

E(5,To) CE(5,T^), if To^T^. 

Let -E'oo('^) '■= Uto>o -^(^' -^o)- Due to the inclusion above we have to check that 

fii{Bl{Mi)\E^{6))=0. 

Denote 

7^(iV) := U {I e Bl{M^) : ■ L = 0}, 

LgZ™o\{0},|L|^Ar 

where = {Wi{I) , . . . ,Wrao{I)) . Let us write (^'"o) as a Fourier series 

= Elsz^o i^fc^e*^-^"" , where F^^ = F^{I""^). Then there exists A/'o > 

such that 



< -, A;=l,---,no. 



Arguing as in the proof of Lemma 2.2, we see that if / ^ TZ{Nq), then 

I y 1 r°RV^-^™°*rft ^-f inf |L-iy-''|r' y 

\ ^ To Jr. ^ Tq\o^\l\^No^ 7 I ^ I 

where W^^ = W^™°(J). The r.h.s of the above inequality can be made smaller than 6/2 
by choosing Tq large enough. So we have 

Bl{M^)\n{No)cE^{6), 

and it remains to show that 

/x,(7^(iVo)) = 0. 
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By Lemma 1.2, 

iy'"°(/)-L^O, VLgZ™»\{0}, 

Since W^°{I) is analytic with respect to / and yU/ is a non-degenerated Gaussian 
measure, then due to Theorem 1.6 in [13], for any L G Z"*", we have 

/i/({/G/i?:Vr"»(/)-L = 0}) = 0. 

Therefore, 

/i,(7^(iVo)) = 0. □ 

Lemma 4.5. Fix any 5 > 0, p > 0. Then for every z/ > we can find Tq > such that 

/i/,^(ro\Ar) < z/, 

whereiV = iV(g^,e,5,To). 

Proo/: Let us denote Te = E{6,To) x T°°, Fi = 5^(Mi) x and Ff := 
(Jjj^^Q F£;(5, To). Since the sets F£;(5, Tq) are increasing with Tq, then from Lemmas 
4.3 and 4.4 we know that 

hm fiijT^ \ F^(5, To)) = fijjT^ \ F^) = 0. (4.16) 

Tq— ^oo 

Let dfii be the measure dfidt on /i^ x M, and 5* ^ be the flow of the perturbed KdV 
equation (II. 2p on h^. We now define following subset of x M: 

B' = [{v,t) : Sliv) G n-;(Fi \ F^(5,To)),t; G i?;(V^),t G [0, i]}. 
By Theorem 3.6, there exists a constant C2(Mi) depending only on Mi such that 

/ii(5') = y /i(5-*(n,-i(Fi \F^(5,To)) f|n,-^(Fo))dt 

< V^(*^^V(n-^(FAr^(5,To))) 

= ^e^^(^^W(FAr^(5,To)). 
For V G n7^(Fo), we define 

S{I,^) = S{v) = {te [0,e-'] : Sliv) G 5;( V^) \ n7;^(F^(5, To))}. 
By the Fubini theorem, we have 

m{B') = [ Mes{S{v))fi{dv), 

Thus 

/i,,^(Fo \ N) = fiiJiil.ip) G Fo : Mes(^(/,y.) > ^r^}) 

^^^^^^^^i^ p,,,(FAr^(5,To)), 

by the Chebyshev inequality. In view of fl4.16p the term on the right hand side becomes 
arbitrary small when To is large enough. The statement of Lemma 4.5 follows. □ 



(4.18) 
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We pass to the slow time r = et. Let v'^ij), r G [0, 1], be a solution of the equation 
1^ and (P(r),<^^(r)) = Ujjv^r)). 

By Lemma 2.1 and (13.31) . we know that for any p ^ 0, the mapping 

where {F){J) = ((Fi)(J), (F2)(J), . . .) is analytic. Hence, there exists C^i^Mi) such that 

\Fj{Ji) - Fj{J2)\p ^ C^{M,)\J, - J2|p, Ji, J2 e Bl{2A'h). (4.17) 

Using Picard's theorem, for any Jq G B^^Mi) there exists a unique solution J{t) of the 
averaged equation (10. 7p with J(0) = Jq. We denote 

T(Jo) :=inf{r>0: |J(r)|p > 2Mi}. 

Now we are in a position to prove the assertion (i) of Theorem 0.2. 

For any p > 0, there exist ni such that 

|F(J,^) - F"H/,^)|p < le-'^^^^^^), (J,^) G 5;(2Mi) x T^, 

|(F)(J) - (F)-(^)lp < ^e-^«(*^^), J G i?;(2M0. 
Choose po such that 

ni 

8$:/^''Po = p-e-^^(^^^ 
i=i 

By Lemmata 4.1 and 4.2, there is a set Fp = A^( g^^° ^ , e, ^j^""^)) < 1/2, such 
that if e is small enough and (/"(O), ^9^(0)) G Fp, then 

4(r) = 4(0)+/ (Ffc)(/(s))ds + 6(r), |efc(r)|<po, rG [0,1], 

for = 1, ■ ■ ■ , ni. Therefore, by ( KIT} and fliTTS]) . 

|/(r)-J(r)|p^ C,{M^)\I{r)~J{T)\,ds + Ur). leo(r)| ^ ^e'^^^^^), 

for (J(0),¥P(0)) G Fp, /(O) = J(0) and |r| ^ min{l, T( J(0))}. By Gronwall's lemma, 

- J{r)\, ^ p, |r| ^ min{l, T( J(0))}. 

Assuming that p << Mi, we get from the definition of T(J(0)) that T(J(0)) is 
bigger than 1. This establishes inequality (10.91) . From Lemma 4.5 we know that 
lim,^o P/,¥p(ro - Tp) =0. □ 

4-2. Proof of the assertion (ii) 

It is not hard to see that the assertion for any ^ Ti < T2 ^ 1 would follow if we can 
prove it for Ti = 0, T2 = 1. So we assume that Ti = 0, and T2 = 1. For any (m, n) G N^, 
we fix a < 1/8, and denote 

Bm{e) := il G B'M) : inf |4| < e"|, 

7^^m,n(e) := U {/ G fii(Mi) : |H^(/) ■ L\ < e"}}. 
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T 

J. r 



Denote 



S{e,m,n,Io,fo) = {r e [0,1] : r(r) G T„,„(e)} 
and fix any z/ > 0. Then using Theorem 3.6 and arguing as in Lemma 4.4 and Lemma 
4.5, we get that, for any (m, n) G N^, there exists open subset FJ^''" C Fq, e^.n > and 
a positive function Pm,n(e), converging to zero as e — )■ 0, such that 

fJ'iA^o - riT'") < ^ and Mes(5'(e, m, ra, Jq, v^o)) ^ Pm,n(e), 
if (/o,^o) eF-'" ande^e„,„. Let 



(m,n)eN2 



then 



(4.19) 
(4.20) 



The sets Ty may be chosen in such a manner that 

Fj,^ C F^2, if z/2 < z/i. 
For any (Jcv^o) G F^, consider a solution (/""(r), v5^(t)) such that 

(r(0),y.^(0)) = (/o,¥.o). 
Fix m G N, take a bounded Lipschitz function g defined on the torus C T°° such 
that Lip{g) ^ 1 and |(7|l^ ^ 1. Let Xlsez™ Fourier series. Then for any 

p > 0, there exists ra, such that if we denote (jn = S|s|^n 5'^^*'"^' 



then 



For any (/o,v9o) G F^, we consider the set S'(e, m, n, /q, v^o)- It is composed of 

open intervals of total length less than T = Pm,n(e)- Proceeding as in Lemma 4.1 and 

Corollary 4.2, we find that for e small enough we have 
"1 



g{^'^^{T))dT 



g{ip)dip 



<p. 



That is 



g{^)fi'f^ f^{d^) - I g{^)d^ 



as e — 7- 0, 



(4.21) 



for any Lipschitz function as above. Hence, p^^^ j,^ converges weakly to dip (see [9]). This 
proves the required assertion with F;^ replaced by F,^. Let us choose 



Then 



v>0 



0, 



by fl4.19|) and fl4.20p . and for any (Jo,(/9o) ^ the required convergence of measures 
holds. This proves the second assertion of Theorem 0.2. □ 
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5. Application to a special case 

In this section we prove Proposition 0.3. Clearly, we only need to prove the statement (ii) 
of assumption A. Let J-" : H"^ ^ M be a smooth functional (for some m ^ 0). If u{t) is 
a solution of ( 10. ip . then 

j^Hu{t)) = (yHnm-V{u) + ef{x)). 
In particular, if J^{u) is an integral of motion for the KdV equation, then we have 

(VJ^(M(t),r(M)) =0, so 

jF{u{t)) = elynumf{x)). 
Since ||m(0)||q is an integral of motion, then 

dt 



^|Kt)||^ = 26(«,/(a;))^e(||«||^+||/(x)||^). 



Thus we have 

|Kt)||^^e^*(|KO)||^ + et||/(x)||^). (5.1) 

The KdV equation has infinitively many integral of motion J'm{u), m ^ 0. The 
integral J7m can be writen as 

m „ 
r=3 m 

where the inner sum is taken over all integer r- vectors m = (mi, . . . ,mj.), such that 
^ rrij ^ m — 1, j = 1, . . . , r and mi + ■ ■ ■ + rrir = 4 + 2m — 2r. Particularly, 
Jo{u) = \ \u\\l. 
Lets consider 

/ = / ■ ■ ■ /("^'^ ■ ■ ■ ^("'■i^dx, mi + ■ ■ ■ + m,, = M, 



where ri ^ 2, M ^ 1, and ^ m^ ^ yU — 1. Then, by Holder's inequality, 

|/|^|M'"^^IIl,,---||/(x)|U,.---|K-IU,,, p, = — ^oo. 

rrij 



Applying next the Gagliardo-Nirenberg and the Young inequalities, we obtain that 

\I\^6\\u\\l + Cs\\u\\^\ W6>0, (5.2) 

where Cs and Ci do not depend on u. Below we denote C a positive constant 
independent of u, not necessary the same in each inequality. Let 

m 

r=3 m 

where mi + ■ ■ ■ + = 6 + 2m — 2r. Using 05. 2p with a suitable 5, we get 

h ^ \\u\t + C\\uf,^ ^ \\u\\l + C{1 + llnlir) + II/IIL (5.3) 
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and 



If u{t) = u{t,x) is a solution of equation flO.ip . then 

It 
1 



I- J'r.(M) = (V J-^(m), e/) ^ e| |m| 1^ + eC(l + I |m| 1^) + e" 



ml 



- C(l + ^ Jn.{u) ^ 2\\u\\i^ + C(l + \\u\ 



2' 

Denote = C(l + ||m(0)||^'") + then from and above, we deduce 



thus 



so 



^( Jm(M) - Cm) ^ ^e(Jm(M) - C„ 



|Kt)||^^4|K0)||^e^^* + C^. 
This prove Proposition 0.3. □ 



Appendix 

Consider the following system of ordinary differential equations: 

X = Y{x), x{0) = xo G M", 

where Y{x) = (Y'i(x), ■ ■ ■ , Yn{x)) : — t- R" is a continuously differentiable map. Let 
F{t,x) be a (local) flow determined by this equation. 

Theorem A (Liouville). Let B{xi, ■ ■ ■ be a continuous differentiable function on 
M". For the Borel measure dfi = B{x)dx in and any bounded open set A C MJ^, we 
have 

where T > is such that F{t,x) is well defined and bounded for any t G {—T,T) and 
X e A. 

For B = const this result is well known. For its proof for a non-constant density B 
see e.g. 
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